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ABSTRACT 

Laporte et al. (2011 1 reported a very high redshift galaxy candidate: a lensed J— band 
dropout (A2667-J1). Jl has a photometric redshift of z = 9.6 — 12, the probability 
density function for which permits no low or intermediate z solution. We here report 
new spectroscopic observations of this galaxy with VLT/XShooter, which show clear 
[Oiii]A5007A, Lya, Ha, and H/3 emission and place the galaxy firmly at z = 2.082. The 
oxygen lines contribute only ~ 25% to the H— band flux, and do not significantly affect 
the dropout selection of Jl. After correcting the broadband fluxes for line emission, we 
identify two roughly equally plausible natures for A2667-J1: either it is young heavily 
reddened starburst, or a maximally old system with a very pronounced 4000A break, 
upon which a minor secondary burst of star formation is superimposed. Fits show that 
to make a 3<r detection of this object in the B— band (V— band), imaging of depth 
AB=30.2 (29.5) would be required - despite the relatively bright NIR magnitude, 
we would need optical data of equivalent depth to the Hubble Ultra Deep Field to 
rule out the mid-z solution on purely photometric grounds. Assuming that this stellar 
population can be scaled to the NIR magnitudes of recent HST/WFC3 IR-selected 
galaxies, we conclude that infeasibly deep optical data AB~ 32 would be required for 
the same level of security. There is a population of galaxies at z w 2 with continuum 
colours alone that mimic those of our z = 7 — 12 candidates. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: individual: 
A2667-J1 — galaxies: star formation 



1 INTRODUCTION 

Broadband imaging in the near infrared (NIR) has enabled 
the detection of candidate galaxies at redshifts between 
z = 7 and 12, using both recent survey instruments on the 



Hubble Space Telescope ( 


McLure et al. 


2010 Finkelstein 


et al. 


2010 |Bouwens et a 


L.||2011b| |Trenti et al.||2011| |Yan 


et al. 


2011 


1 and wider-but-shallower observations from the 



al in prep). With specific astrophysics in mind, these cata- 
logues may be used to infer the cosmic history of star forma- 
tion, the photon budget available for reionization, and also 
as an independent test of the re-ionization edge through the 



* Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 087.A-0118 
f E-mail: matthew.hayes@irap.omp.eu (MJH) 



prevalence of the Hi Lya emission line ( Pentericci et al. 2011 
|Schenker et al.|2012| [Curtis-Lake et al.||2012j 

It is vital, therefore, to understand the lower-z objects 
that could contaminate such samples, and the frequency 
with which they do. At the highest z, contamination rates 
have been estimated using population synthesis methods 
( Trenti et al.|2011 Bouwens et al.|2011b |, and shown to be 
non-negligible. Empirical studies at z ~ 2 have also found 
a population of galaxies with emission lines strong enough 
to mimic a strong spectral break, suggesting roughly similar 
contamination levels ( |Atek et al.||20llj ). Yet the empirical 
interloper rate, measured directly by deep spectroscopy in 
the z > 7 samples, remains broadly unknown. 

To these ends our group has been carrying out its own 
multi-wavelength studies, taking advantage of gravitational 
lensing to search for intrinsically fainter background objects 
(Laporte et al. 2011 Lll) Originally, Lll selected ten lensed 
galaxies that satisfy criteria for z—, Y-, and J— dropout 
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galaxies. With deep optical, near-infrared and Spitzer/IRAC 
photometry, traditional photometric redshift (phot-z) meth- 
ods placed all galaxies in the range z = 7 — 12. However two 
of the Lll sources were found to he at z ~ 2 ( |Boone et al 



20111, hinting that strong optical lines may have mimicked 
a Lyman break. In cases where very deep optical imaging 
is available, this source of contamination is thought to be 
small, as the line strengths imply lowish masses/metallicities 
that in turn should exhibit blue continuum colours and re- 
main detectable, particularly when stacked. However other 
sources of contamination exist, specifically relating to dust 
attenuation and strong 4000A breaks: in both of these cases 
an intrinsically red spectrum in the restframe UV and opti- 
cal may cause a galaxy to go undetected in all bandpasses 
bluer than the first dropout band. 

In this Letter we present the optical and NIR spec- 
trum of a bright J— band dropout galaxy (A2667-J1) - phot- 
z = 11. 9^2 3' 1 , no low- z solution and no Herschel detection 
- obtained with VLT/XShooter. In Lll we were skeptical 
about the very high-z interpretation because of the large 
implied luminosity, and the detection of a faint source in 
HST/zgso imaging (although the z — H colour still exceeds 
2 magnitudes). However the spectrum shows several strong 
emission lines that are consistent only with z = 2.082, im- 
plying it is an interloper. In Sect. [2] we present the new 
observations and measurements of emission lines. In Sect. [3] 
we discuss inferences we can draw from the emission lines, 
and explanations for the strong spectral break and selection 
function, and re-model the stellar continuum. In Sect. [4] we 
discuss the implications that this galaxy may have upon 
other very high-z surveys. All magnitudes are AB. 



2 OBSERVATIONS AND MEASUREMENTS 
2.1 X-Shooter observations 



A2667-J1 (q=23:51:34.9, 5= 



with X-Shooter (D'Odorico et al. 2006 



26:3:32.7, J2000) was observed 
Vernet et al.|2011[ ), 



mounted at UT2 of ESO/VLT, on the night of 24 August 
2011. The total integration time was short, but at 4 x 900 s 
was sufficient to find several strong emission lines. We used 
the 1" (0.9") slit for the UVB channel (VIS and NIR), pro- 
viding an almost constant aperture between the U— and 
K— bands. More importantly the NIR spectral resolution of 
R > 5000 is sufficient to resolve the atmospheric OH airglow 
lines. Between each observation, the object was dithered 
along the slit by 2.5 ". The average atmospheric seeing was 
0.76 ". Data were reduced and combined using the EsoRex 
X-Shooter pipeline version 1.3.7, and one dimensional spec- 
tra were extracted from the frames using the APALL task in 
noao/iraf with an aperture width of 1.5 ". 



2.2 Emission lines in the spectrum 

We clearly identified a strong emission line at wavelength 
A = 15437A (Fig. [j] - after masking regions contaminated 
by a nearby airglow line and fitting a Gaussian profile, 
we measure a flux of (18.8 ± 0.96) x 10" 18 erg s" 1 cm -2 . 
The //-band magnitude is 25.21 ± 0.08 (Lll), and as- 
suming the quoted bandpass parameters (A c = 1.620um, 
FWHM=2890A), we compute an observer frame equivalent 



Table 1. 


A2667-J1 fluxes in 


CLLLlBoHJII 




opccics 




^mcas 


z 


r 1UX 




A 


A 




10 ^ s ^ cm ^ 


Hi Lya 


1215.67 


3743.0 


2.079 


9.41 ± 2.25 


Hen 


1640.42 






< 4.87 


[On] 


3727.09 






< 2.72 


[On] 


3728.79 






< 2.72 


Hi H/3 


4862.72 


14998.2 


2.084 


7.77 ± 4.22 


[Oiii] 


4960.30 


15288.2 


2.082 


7.06 ±2.72 


[Om] 


5008.24 


15436.6 


2.082 


18.8 ±0.96 


Hi Ha 


6564.61 


20234.6 


2.082 


22.7 ± 4.08 


[Nil] 


6583.46 


20290.2 




< 3.62 



width, EW, of (679 ± 71) A. The wavelength of the observed 
emission line falls close to the peak transmission so no filter 
wing corrections are necessary. 

This high EW restricts the possible line identifica- 
tions - after testing various possibilities, we determine it 
to be [Olll]A5007A and also identify Lya, H/3, and Ha. This 
places the redshift of Jl firmly at z = 2.082. After the re- 
moval of a nearby OH line [Oni]A4959A is also visible but 
no HeiiA1640A, [Nn]AA6548,6584, or [On]AA3726,3729A are 
visible in either the ID or 2D spectra. Example spectra are 
shown in Fig. [I] and the measured line fluxes (or limits) are 
listed in TableTl] 



3 WHAT IS A2667-J1? 

3.1 Basic inferences from the nebular lines 

We first discuss a few line ratios. Both Ha and H/3 can be 
measured, although the significance of H/3 is low (2a). The 
Ha/H/3 ratio is 2.92 ± 1.67, consistent with the unobscured 
value of 2.86 (T = 10 4 K), which corresponds to Eb-v = 
± 0.5 (la, Calzetti et al.|2000 l. Assuming no obscuration, 
the star for mation rate (SFR) inferred from Ha is 4.1 ± 



0.7M o yr ( Kennicutt 



1998 



corrected for a magnification) , 
which is towards the lower end of the envelope of [Hayes et al. 
(2010b I for z ~ 2 Ha-selected galaxies. At the la upper 



limit of dust attenuation the intrinsic SFR will be a factor 
of 5 higher. Regarding Lya, assuming pure recombination 
under case B and the permitted range of Eb-v, Lya escape 
fractions falls in the range 1-5%, again close to the z ~ 2 
average ( Hayes et al.|201Qa i. 

Neither side of the [On] doublet is detected (< 
1.92 erg s _1 cm -2 ), despite the region being free of both tel- 
luric and airglow features. This allows us to put a reasonable 
limit on the excitation parameter, P [=[Om]/([Oni] + [Oll])] 
of 1-OIq j_, which is not greatly affected by the range of per- 
mitted dust contents. P is known to correlate strongly with 
the ionization parameter, and the hardness of the ionizing 
spectrum (therefore the temperature of the stellar popula- 
tion, E.g. Pilyugin 2001 1, and the spectrum suggests that the 



stellar population of A2667-J1 is still hot enough to keep all 
the oxygen doubly ionized. We may also calculate the /?23 
index [([Oll]±[Om])/H/3] which is not particularly high at 
3.2±1.6; the high value of P, the low R23, and the absence of 
a clear [Nil] lines surrounding Ha suggest that we are look- 
ing at a relatively low metallicity burst of star formation. 
From the N2 index, metallicity 12 + log(0/H) < 8.39 (la, 
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Figure 1. Examples of the extracted ID spectra around strong detections. Left shows Hi Lyct, centre shows [Oin] A5007A, and right 
shows Hi Ho. From top to bottom each panel shows the 2D spectral images, the ID spectral extractions, the sky transmissivity /telluric 
absorption, and the OH airglow spectrum. Lyo presents with its characteristic asymmetric shape. [Om] is the strongest line, and detected 
with S/N = 20, while Ha appears somewhat contaminated by both telluric bands and skylines. See text for details. 



|Pettini fe~P agcl 2004]), which would place Jl roughly mid- 
way through the distribution of |Erb et al.| ( 2006| ) at z ~ 2. 
While Jl is very compact, the [Om]/H/3 ratio and [Nil] /Ha 
limit are fully consistent with a nebula produced by star- 
formation, not an active nucleus (Kauffma nn et aL||2003[ ). 



3.2 Line contribution and colour selection 

A2667-J1 met the selection criteria for a z ~ 11 galaxy in 
[(J — H) vs. (H — K s )] colour space denned in Lll and op- 
timized for the removal of stars and lower z galaxies. It is 
important to note that every survey has its own combina- 
tion of wavelength coverage and depth, and therefore the 
selection criteria will naturally differ. Indeed most compa- 
rable studies are based upon HST data, and consequently 
include no K— band criteria for their J— band dropouts. De- 
tails aside, it is well known that [Om] and Ha lines may 
have EW high enough to dominate broadband fluxes (E.g. 
Schaerer fc de Barros|2009l |Trump et al.|2011 1 causing mid- 
z galaxies to interlope into very high-2 samples (Atek et al 



2011 1. A vital question, therefore, is whether Jl was selected 
because of this. Combined [Olll]+H/3 provide EW w 900 A, 
and these lines contribute just 25% to the //—band flux 
(FWHM = 2890 A). If we subtract this nebular contribu- 
tion, the colour of the (J — H) spectral break changes from 
the measured value of > 2.3 magnitudes to > 2.0, and still 
easily clears the dropout criterion of (J — H) > 0.76. After 
subtracting Ha from the K s band we recompute the (H—K s ) 
colour as 0.50, only slightly redder than the 0.44 measured in 
Lll and still on the border of the selection box. Therefore 
even without the strong nebular lines we would still have 
selected Jl as a very high-z candidate. 

There is the need to explain a spectral discontinuity 
with amplitude > 6 in /„ between the J— and //—bands 
(restframe 5250 and 4080 A), and amplitude > 10 between 
Y (restframe 3310 A) and H. This could result from either 
heavy dust obscuration or a strong 4000A break - the for- 



mer seems inconsistent with Hq/H/3, while the latter seems 
inconsistent with even the presence of strong emission lines. 



3.3 (Re-) Modeling the stellar population 

Before modeling, we take advantage of new Spitzer data re- 
duction methods and re-perform the IRAC photometry us- 
ing the latest pipeline processed data. The high-level sci- 
ence data from the Lensing Cluster GTO Survey (Program 
83, PI. Rieke) were re-obtained, and photometry performed 
in 1 arcsec circular apertures paying special attention to 
the sky properties (level, noise, skewness) in nearby empty 
apertures. In Table [5] we present updated stellar (i.e. line- 
subtracted) VLT photometry, updated IRAC photometry, 
and the original HST photometry from Lll. We note that in 
the optical and NIR colours, A2667-J1 appears rather sim- 



ilar to the unusual infrared object HDFN-JD1 (Dickinson 
et al. 1120001), although it is 0.7 magnitudes bluer in H — K. 



- the nature of that source also remains unclear. 

With new stellar (i.e. line-corrected) photometry, we 
model the stellar population with Hyper-z ([Bolzonclla et al. 
2000[), using the CB2007 modification of the |Bruzual fc 



Chariot ( 2003[l population synthesis models to generate evo- 



lutionary templates. We adopt a Salpeter IMF, lock the red- 
shift z — 2.082, and fit dust attenuation, star formation his- 
tory (SFH), and time since the onset of star formation. The 
range of included SFHs include exponentially decaying SFRs 
with e-folding timescales as listed in the abscissa of Fig [2] 
panel (b). For comparative purposes we perform fits at three 
different metallicities: Z = 0.0001, 0.004, and 0.02. We show 
the results in the left panel of Fig. [2] with histograms of 10 4 
Monte Carlo realizations in the right panel. 

Firstly we note that for the three different metallicities, 
the values of \ 2 that describe the best fit are effectively con- 
sistent, and we do not have the power to distinguish between 
them. Metallicity does however modulate the frequency with 
which other results are obtained, and there are two differ- 
ent and entirely degenerate ways of reproducing the break: 
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Figure 2. SED fitting results using Hyper-z and the CB2007 synthesis models. Panel a. (left) shows the fits, with data points and 
limits as solid black boxes; error bars are la and horizontal bars represent the FWHM of the filters. Open boxes show the IRAC points 
when corrected for Paschen and Brackett line emission. The solid coloured lines show different stellar metallicities; inferred quantities are 
displayed in the corresponding colour. The right panel shows the distribution of the quantities recovered from fits to 10 4 Monte Carlo 
realization of the data. The same metallicity colour scheme applies. We show SFH r (b.), age (c), Eb—v (d.), and stellar mass (e.). 



Table 2. Updated stellar fluxes, corrected for emission lines 



are rejected by three arguments: (i.) the necessary attenua- 



Telescope 


obs A c 


rest A c 


magnitude 


& bandpass 


[ l"n ] 


[ H m ] 


[AB] 


VLT/FORS2/I 


0.791 


0.257 


27.39 ±0.49^1 


HST/ACS/F850LP 


0.911 


0.296 


27.39 ±0.18""^ 


VLT/FORS2/z 


0.914 


0.297 


> 27.72 a 


VLT/HA WKI/Y 


1.03 


0.334 


> 28. l a 


VLT/HAWKI/J 


1.26 


0.409 


> 27.5 a 


VL T/HA WKI/H 


1.62 


0.526 


25.53 ±0. lip] 


VLT/HAWKI/K S 


2.15 


0.698 


25.02 ±0.16^ 


Spitzer/IRAC/Chl 


3.58 


1.16 


24.48 iO.lfl 


Spitzer/IRAC/Ch2 


4.53 


1.47 


24.70 ± 0.24^ 


Spitzer/IRAC/Ch3 


5.78 


1.88 


23.65 ±0.41 c 


Spitzer/IRAC/Ch4 


8.05 


2.61 


> 24.62 c 



as reported in Lll 
emission line subtracted 
photometry modified from Lll 



either with a large dust content or an old population with 
a strong 4000A break. The latter scenario would imply a 
dual stellar population: one evolved enough for a strong 
break, the other young enough to produce the emission lines. 
As usual a strong 4000A break is more frequently seen at 
higher metallicities, where stellar evolution runs faster. In- 
deed models at extremely low Z prefer the young dusty so- 
lution in > 80% of cases, models at Z = 0.004 seem to be 
rather evenly split between the young/dusty vs. old solu- 
tions, and the Z — 0.02 models almost exclusively prefer 
ages above 2 Gyr and little dust. The Z = 0.0001 models 



tion is higher than permitted from Ha/H/3 (Sect 3.1 1, while 
usually nebular Eb-v exceeds that of the stars (Calzetti 
|et al.|[2000| |; (ii.) such metallicities should produce higher 
[Oiii]/H/3; and (Hi.) a dusty metal-poor galaxy seems an un- 
likely combination. However the Z — 0.004 models all seem 
permissible, and we do not have the information to formally 
distinguish between the two scenarios. 

We note that in the case of the composite population, 
the recent burst must also have a stellar population that has 
not been subtracted from the photometry prior to fitting. In- 
deed its contribution to the J— band flux is not permitted 
to be large, otherwise the J — H break would be dimin- 
ished. To investigate this we take the restframe Ha EW 
(280A) and compute the range of ages at which this EW is 
exceeded: the star-forming population must have an age be- 
low w 5Myr (|Leitherer et al.||1999] Starburst99). Over this 



age range we then compute the range of permitted J mag- 
nitudes and without reddening the population will remain 
completely undetected in J for ages below 2 Myr. However 
the Ha EW is independent of reddening while the J — K s 
colour is not; the la limit of Eb-v allows up to a factor of 
3 greater attenuation at J than at K and ages up to 4 Myr 
still predict a J— band non-detection. On the other hand, a 
rather blue continuum may explain the HST 2sso detection. 
For this range in parameter space, the burst contributes be- 
tween 0.4 and 4 % of the total stellar mass. 

Although the IRAC channels show the largest errors, 
the 5.8um channel carries a residual of « 1.2cr. We note 
that the higher metallicity stars show the strongest 1.6u.m 
stellar feature, which is probably the cause of the marginal 
preference for these models. However the IRAC points also 
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contain the Paschen and Brackett lines, for which we also 
make a first order attempt to correct (assuming f0 4 K gas). 
Since the nebular extinction is not well constrained, we make 
two corrections, one derived for Eb-v=0 and the other for 
Eb-v=0.5. These modified points are shown in Fig [2] and 
the largest correction (where Paa falls in 5.8(J.m) amounts 
to just 0.06 magnitudes - consequently the corrections in 
channels 1 and 2 are not visible. We do not re-perform the 
SED fit, but note that this correction does just bring the 
channel 3 data point to about la of the fit. 

Finally we examine the phot-z solutions that are al- 
lowed by the new line-free photometry, and we repeat the 
fits with the redshift free. The peak at z = 10 — 12 is still 
recovered but now represents the tertiary solution. The pri- 
mary solution lies at z = 2.6 — 3.3 which is actually in good 
agreement with the secondary peak from Lll (z = 2.8), and 
the secondary grouping prefers z < 1.5. However, more inter- 
estingly, only 0.2% of the realizations have phot-z within la 
of the spectroscopically measured value, and it is not easy 
to obtain a good SED fit at the measured redshift when 
the templates comprise only single stellar populations. This 
lends support to the scenario of multiple populations, and 
the hypothesis that the zsso and NIR detections do not see 
the same stars. 



4 IMPLICATIONS FOR HIGH-Z SURVEYS 

A2667-J1 can be described in two ways, but regardless of 
the true population, one salient fact remains: was photo- 
metrically selected as a z ~ 11 candidate, but actually lies 
at z w 2. 

Because of our homogeneous coverage across the NIR, 
Jl was selected in two colour space: (J — H) and (H — K s ). 
The selection, therefore, is more similar to lower z meth- 
ods (e.g. z— and Y— dropouts) and somewhat different from 
other J— drop searches (Henry et al. 2007 Bouwens et al. 



2011a), which lack comparably deep if— band data and 



therefore employ a one-colour selection. The requirement of 
(H — K s ) < 0.5 (a blue UV slope) is specifically intended to 
fortify the selection against the inclusion of mid-z galaxies 
with strong optical breaks, but a z « 2 galaxy is precisely 
what we have found. On the other hand, a non-detection in 
bandpasses bluewards of the break is also required by typical 
surveys to reject lower- z galaxies and, while Jl is undetected 
in ground-based z, Y, and J imaging, it is detected in HST 
Z850 (albeit with a colour similar to the dropout criterion). 
For Jl, imaging bluer than the z band would be effectively 
useless in evaluating its status as an interloper: from the 
best-fitting SEDs and Monte Carlo simulations described in 
Sect. [373J we calculate that in the HST F435W and F606W 
(B and V) we would need to reach depths of AB=30.2 and 
29.5 in order to detect the continuum. In other words, this 
relatively bright //—detected J— drop would need imaging 



with the depth of the Hubble Ultra Deep Field (Beckwith 



et al. 2006 1 to be rejected as very high-z. A magnitude 



H = 28 galaxy of these colours in the HUDF or CANDELS 
surveys would still remain undetected by in the available 
HST z— band imaging, and require imaging reaching V = 32 
in order to detect the bluer continuum. It currently remains 
unclear how abundant galaxies of this type are, and thus 



systematic campaigns of spectroscopic followup will remain 
a vital pursuit, even at the faintest magnitudes. 
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